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ABSTRACT 

We construct a protostellar disk model that takes into account the combined effect of viscous 
evolution, photoevaporation and the differential radial motion of dust grains and gas. For T Tauri 
disks, the lifetimes of dust disks that are mainly composed of millimeter sized grains are always 
shorter than the gas disks’ lifetimes, and become similar only when the grains are fluffy (density 
< 0.1 g cm -3 ). If grain growth during the classical T Tauri phase produces plenty of millimeter 
sized grains, such grains completely accrete onto the star in 10 7 yr, before photoevaporation 
begins to drain the inner gas disk and the star evolves to the weak line T Tauri phase. In the 
weak line phase, only dust-poor gas disks remain at large radii (> 10 AU), without strong signs 
of gas accretion nor of millimeter thermal emission from the dust. For Herbig AeBe stars, the 
strong photoevaporation clears the inner disks in 10 * * 6 yr, before the dust grains in the outer disk 
migrate to the inner region. In this case, the grains left behind in the outer gas disk accumulate 
at the disk inner edge (at 10 — 100 AU from the star). The dust grains remain there even after 
the entire gas disk has been photoevaporated, and form a gas-poor dust ring similar to that 
observed around HR 4796A. Hence, depending on the strength of the stellar ionizing flux, our 
model predicts opposite types of products around young stars. For low mass stars with a low 
photoevaporation rate, dust-poor gas disks with an inner hole would form, whereas for high mass 
stars with a high photoevaporation rate, gas-poor dust rings would form. This prediction should 
be examined by observations of gas and dust around weak line T Tauri stars and evolved Herbig 
AeBe stars. 

Subject headings: accretion, accretion disks — circumstellar matter — planetary systems: formation — 
solar system: formation 


1. Introduction 

Theoretical models for the evolution of proto¬ 
stellar disks have to satisfy a number of observa¬ 
tional constraints regarding the relative timescales 
on which their dust and gas are dispersed. These 
constraints are: 

(a) The majority of young stars of spectral type 
GKM exhibiting evidence of gas accretion onto the 
stars (i.e., the classical T Tauri stars, henceforth 
CTTSs; age ~ 10 6 yr) are observed to possess dust 
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disks of ~ 100 AU in size as evidenced by their mil¬ 
limeter emission. For example, Osterloh & Beck¬ 
with (1995) detected 32 of the 62 CTTSs (52%) 
in their sample at a typical 3er level of 10 — 20 
mJy (see also Beckwith et al. 1990). Thus the 
dust in the outer disks cannot disperse below the 
present detection limit on timescales much shorter 
than that of the gas in the inner disks. 

( b ) By contrast, the majority of stars of spectral 
type GKM which lack spectroscopic accretion in¬ 
dicators (i.e., the weak line T Tauri stars, WTTSs) 

are usually not detected at millimeter wavelengths 

[5 of 46 WTTSs (11%) in the sample of Osterloh 

& Beckwith 1995; 2 of 10 WTTS in Dutrey et 
al. 1996], even in recent high sensitivity surveys 
achieving 3<r upper limits of a few mJy (1 of 12 in 

Duvert et al. 1999). Thus the dust in the outer 
disks cannot disperse on timescales much greater 
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than that of the gas in the inner disks. 

At first sight, the simplest explanation of (a) 
and ( b ) is to invoke a constant dust-to-gas ratio 
both in space and in time, as is usually assumed 
in modeling the disk evolution over 10 6 yr (e.g. 
Hartmann et al. 1998; Armitage et al. 1999, 2003; 
Clarke et al. 2001). If most grains are as small 
as ~ 1 n m, the dust and the gas are well coupled 
and so the constant dust-to-gas ratio over 10 6 yr is 
probably a reasonable assumption. However, such 
an assumption is questionable for dust grains of 
millimeter size or above, given their radial migra¬ 
tion (Takeuchi & Artymowicz 2001; Takeuchi & 
Lin 2002, 2003, and 2005, henceforth TL02, TL03, 
and TL05, respectively). Figure 1 of TL05 shows 
that the orbits of bodies in a size range of 1 mm to 
10 nr decay in less than 10 6 yr, which is the typi¬ 
cal age of CTTSs (Kenyon & Hartmann 1995). In 
addition, the constant dust-to-gas ratio does not 
appear to be compatible with a third constraint, 
i.e., 

(c) The transition timescale, over which stars 
lose circumstellar diagnostics, is an order of mag¬ 
nitude less than the lifetime of the disk. This con¬ 
clusion is based on the relative paucity of ‘tran¬ 
sition objects’ with infrared colors characteristic 
of optically thin inner disks (Strom et al. 1989; 
Skrutskie et al. 1990; Kenyon & Hartmann 1995). 
The rapid turn-off of the inner disk, and of gas 
accretion diagnostics, can be reproduced by mod¬ 
els that combine viscous evolution with photoe¬ 
vaporation by the central star (Clarke et al. 2001; 
Armitage et al. 2003). These models however pre¬ 
dict a long lived outer gas disk (at radii > 10 AU). 
If the dust-to-gas ratio in the outer disk remains 
constant over the typical age of CTTSs (~ 10 6 yr), 
such evolution probably violates constraint ( b ), as 
discussed by Clarke et al. (2001). 

In TL05, we showed that the long lived millime¬ 
ter flux of CTTS disks indicates that a sufficient 
amount of small grains less than ~ 1 mm must 
remain in the disks of ~ 100 AU scale for more 
than 10 6 yr. Otherwise, grains larger than 1 mm 
(and smaller than 10 m) rapidly spiral inward to 
the star due to strong gas drag (see Fig. 1 of 
TL05). The population of the dust grains larger 
than 1 mm disappears quickly from the disk, and 
thus the millimeter flux of the dust attenuates be¬ 
low the observational detection limit within 10 6 
yr. On the other hand, the evolution of the spec¬ 


tral index at millimeter wavelengths suggests grain 
growth to at least millimeter sizes (Beckwith & 
Sargent 1991; Mannings & Emerson 1994; Beck¬ 
with et al. 2000; Calvet et al. 2002; Testi et al. 
2003, however see also Agladze et al. 1996; Men- 
nella et al. 1998; Dupac et al. 2003 for other 
interpretations). Combining these two considera¬ 
tions, we concluded that only a narrow range of 
the grain size around 1 mm is responsible for the 
observed millimeter fluxes of CTTSs. 

Based on the above result, we further inves¬ 
tigate the relation of the lifetimes between the 
outer dust and inner gas disks. In this paper, we 
show that, using plausible assumptions about the 
grain population, the predictions of viscous pho¬ 
toevaporation models are compatible with all the 
constraints (a)-(c), once one accounts for the dif¬ 
ferential radial motion of gas and dust predicted 
by simple grain drag models. We describe such 
a ‘successful’ model in which we assume (i) that 
the dominant opacity at millimeter wavelengths 
derives from grains of sizes close to 1 mm and (ii) 
that there are no significant sources or sinks of 
grains in that size range (i.e., that one may ne¬ 
glect the net creation or destruction of such grains 
due to coagulation or the collisional break up of 
grains). We set out further observational predic¬ 
tions of our ‘successful’ model. 


2. A ‘Successful’ Model for the Evolution 
of Dusty Disks 


2.1. Evolution of the Gas 


We model the evolution of a gas disk due to 
the combined action of viscosity and photoevap¬ 
oration by ionizing radiation from the central 
star. We follow the models of TL05 and Clarke 
et al. (2001). The initial gas density profile is 
Eg = 3.5 x 10 2 (r/AU) _1 g cm -2 with the outer 
radius of r out = 100 AU. The disk has a gas mass 
M g = 2.5 x 10~ 2 M a inside 100 AU. The gas den¬ 
sity obeys the equation, 


-E — 

dt 9 r dr 




— — E„ 


(1) 


The turbulent viscosity is expressed as v t = 
ac 2 /Hk oc r, where c is the sound speed and fix is 
the Keplerian angular velocity. The temperature 
profile, T = 278 (L/Lg^/^r/AU) -1 / 2 K, is as¬ 
sumed to be unchanged during the disk evolution. 
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The sink term E w in the right hand side is the 
wind mass loss rate by photoevaporation, which 
has non-zero values outside the gravitational ra¬ 
dius r g = GM/c 2 ~ 10 14 cm, where M = 1 M 0 
is the central star’s mass and Cj = 10 6 cm s _1 is 
the sound speed of the photo-ionized gas (Shu et 
al. 1993; Hollenbach et al. 1994). The wind mass 
loss rate per unit area is 

t w = 2an 0 nm ■ f , ( 2 ) 


where ?ro is the number density of hydrogen ions 
at the base of the photoevaporating flow, tor is 
the mass of a hydrogen atom. The base density is 

n 0 = 5.7 x 10 4 $l{ 2 r~i/ 2 cm -3 , (3) 

where $44 is the ionizing photon luminosity of the 
star in units of 10 41 s _1 , and r g 44 is the gravita¬ 
tional radius in units of 10 14 cm. We introduce a 
function / to express a drop in the wind mass loss 
rate inside r g , 


f 


1 + exp 



( 4 ) 


where the transition width is assumed as A r„ = 
0.1 r g . 

In our model, the ionizing photon (EUV) flux 
$41 is constant with time. Because absorption 
by interstellar hydrogen atoms prevents us observ¬ 
ing EUV fluxes of young stars, there is no direct 
observational constraint on it, and the origin of 
emitting EUV photons is also unclear. If EUV 
photons are emitted from accretion shock at the 
stellar surface, its flux decays with time and be¬ 
comes inefficient as the disk accretion ceases (Mat¬ 
suyama et al. 2003; Ruden 2004). Alexander et 
al. (2004) also showed that accretion shock can¬ 
not emit enough EUV photons due to absorption 
by the stellar atmosphere and by the accretion 
column. Thus, chromospheric activity of young 
stars is probably responsible for their EUV flux. 
Recent analysis of the far-ultraviolet (FUV) line 
ratios in T Tauri stars (Alexander et al. 2005) 
confirmed ionizing photon production rate in the 
range 10 41 — 10 43 s _1 . In this paper, we simply 
assume a constant EUV photon luminosity $ 44 . 


2.2. Evolution of the Dust 

The initial dust-to-gas ratio, /dust> is assumed 
to be constant with r and equal to 10 -2 . Thus, 


the initial dust density profile is E^ = 3.5 x 
10 2 /dust(r/AU) _1 g cm -2 . This is equivalent to 
assuming that the dust and gas mixture arrives in 
the disk plane in the proportions in which they 
are found in the interstellar medium. Naturally, 
we have no strong evidence that this is indeed the 
case and that some fractionation of the dust and 
gas has not already occurred during the collapse 
of the prestellar core and during the initial evo¬ 
lution of the disk. This issue however is beyond 
the scope of this paper. We instead explore how 
the differential radial motion of dust and gas al¬ 
ters the assumed initial constant dust-to-gas ratio. 
The evolution of the dust surface density subject 
to turbulent diffusion and gas drag is modeled ac¬ 
cording to TL05 as 

v^Erf + [r(Fdif + E d Vd)\ = 0 , (5) 

where F^it is the diffusive mass flux caused by 
gas turbulence, and Vd is the dust radial velocity. 
We adopt the assumption that the diffusive mass 
flux is proportional to the concentration gradient, 
and that the diffusion coefficient is the same as 
the gas viscosity, which is applicable to small pas¬ 
sive particles tightly coupled to the gas, i.e., the 
stopping time of the grains due to gas drag, t s . is 
much smaller than the orbital time, U^ 1 , and the 
Schmidt number Sc = 1. In this case, 


(|) ' (6) 


The radial velocity of the dust subject to gas drag 
is calculated according to equation (23) of TL02 
as 


Vd 


T s 1 v g - T)V K 
T s + TV 1 


Vg - r]T s VK , 


( 7 ) 


where Vk is the Keplerian velocity, T s = UQk is 
the non-dimensional stopping time of the grain, 
is the divergence factor of the gas rotation velocity 
from the Keplerian value. The second equality is 
valid for T s <C 1. The gas radial velocity v g is (eq. 
[11] of TL05) 


r 4 / 2 E ff dr 


(r 1 / 2 vt E g ) . 


( 8 ) 


The non-dimensional stopping time of the grain 
T s that is normalized by the local orbital time is 
(eq. [21] of TL05) 


_ TrppS 

S ~ 2S g 
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where p v is the grains’ physical density and s is 
the grain size. The ratio of the radial gas pressure 
gradient to the gravity is (eq. [20] of TL05 with 
<7 = - 1 / 2 ) 

_ = _L _OH /M 2 (J_^g_ _ 7\ 

r^KPg dr \ r ) V s s dr 4 ) ' 

( 10 ) 

In the above formulation we use midplane val¬ 
ues and neglect the fact that the gas rotation ve¬ 
locity and r) are functions of the height z from 
the midplane and that the stopping time T s also 
varies due to the gas density change with z. This 
treatment is justified if grains sediment to the 
midplane quickly enough. From equation (14) of 
TL03, the timescale of sedimentation is r se a = 
\z/v z ,d\ = (TsOk) - 1 , where v Z)C l is the verti¬ 
cal velocity of the grain. This is smaller than 
the timescale of the grain’s radial drift, Tdust = 

| r/(v d - v g )\ = 7? _1 (T s f2 K ) _1 , by a factor ry < 1 . 
Sedimentation also proceeds more quickly than ra¬ 
dial diffusion, whose timescale is Tdif, r = r 2 /v t ss 
a _ 1 77 _ 1 T s (T s flK) _1 , if a _ 1 ? 7 _ 1 T s > 1. In our 
standard model, a = 10~ 3 , 77 ~ 10~ 2 , and T s > 
10 ” 4 for millimeter sized grains at r > 10 AU 
(see Fig. 2 of TL05), and thus the condition 
a~ 1 rj~ 1 T a 1 is satisfied. Turbulent diffusion of 
moderate strength (a ~ 10 -3 ) cannot lift millime¬ 
ter grains up beyond the disk scale height against 
their sedimentation (for r > 10 AU, see Fig. 4 of 
TL02). 

2.3. Opacity and Grain Size 

In this ‘successful’ model, we consider a popu¬ 
lation of grains with size 1 mm. The justification 
of this choice is that for a grain size distribution 
n(s) oc s~ s with 5 < 4, the opacity at wavelength 
A ~ 1 mm is dominated by the largest grains, pro¬ 
vided that they are small enough to be described 
by the Rayleigh scattering theory (i.e., s < A; see 
TL05). Thus, if the grain size distribution extends 
up to ~ 1 mm in size, the grains can be modeled 
as a single sized population of the largest grains as 
regards their opacity at millimeter wavelengths. 

Evidently, this will not be a good approxima¬ 
tion if the maximum of the grain size distribution, 
s ma „ is either substantially less than or greater 
than 1 mm. We note that the dust opacity index, 
/3, (which measures the frequency dependence of 
the dust emissivity according to k v oc is < 1 


for CTTSs (Beckwith & Sargent 1991; Mannings 
& Emerson 1994; Kitamura et al. 2002), which is 
consistent with grain growth to at least 1 mm. On 
the other hand, if the dust were instead to reside 
mainly in much larger grains, then the opacity at 
1 mm would be considerably reduced (Fig. 4 in 
Miyake & Nakagawa 1993) and could not explain 
the observed millimeter fluxes from CTTSs. 

Finally, in order to compute the disk emis¬ 
sion at 1 mm, we follow Hartmann et al. (1998) 
and Clarke et al. (2001) in assuming that the 
grain opacity is given by k v = 0.1(^/10 12 Hz) 
and that the dust temperature is a fixed power 
law function of the distance from the central star 
[T = 278 (. L/Lq ) 4 / 4 (r/l AU )" 1 / 2 K], 

3. Gas and Dust Disks’ Lifetimes 
3.1. Model T Tauri Stars 

Following TL05 and Clarke et al. (2001), we 
take the parameters as M = 1 M 0 , L = 1 L 0 , 
Tout = 100 AU, M g = 2.5 x 10 - 2 Mq, / dust = 
10 ~ 2 , s = 1 mm, p p = 0.1 — 1 g cnr -3 , a = 10 ~ 3 , 
r g = 10 14 cm, $41 = 1 , and c* = 10 6 cm s _1 , for 
T Tauri disks. The equations for the gas and dust 
disks are solved in the same manner as TL05. 

The dashed lines of Figure 1 depict the evolu¬ 
tion of the gas surface density and illustrate that 
viscous evolution dominates the emptying of the 
disk until the point is reached at which the accre¬ 
tion rate becomes compatible with the mass loss 
rate due to photoevaporation (~ 10 _1 ° Mq yr _1 
for an assumed ionizing luminosity of $ = 10 41 
s _1 ). This mass loss, which is restricted to regions 
of the disk at radii > 7 AU (where photoionized 
material becomes unbound), is eventually able to 
remove all of the mass flux from the outer disk 
and thus cuts off the inner disk from re-supply. 
Consequently, the inner disk empties on the vis¬ 
cous timescale of the disk at ~ 7 AU (~ 10 6 yr for 
this model), which is much less than the inner gas 
disk’s life time of ~ 10' yr. Gas however remains 
at radii > 7 AU for a further 10 7 yr, this being the 
timescale on which gas can viscously diffuse from 
the outermost disk (~ 100 AU) and reach radii 
of ~ 10 AU at which photoevaporation becomes 
effective. 

The solid lines of Figure 1 depict the corre¬ 
sponding evolution of the dust for two cases: (a) 
the grains’ physical density, p p = 1 g cnr ~ 3 and ( b ) 
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Pp = 0.1 g cm -3 . Evidently, from equation (9), 
the lower density ‘fluffy’ grains are more tightly 
coupled to the gas. Since for both grain densi¬ 
ties, T s C 1 (apart from where the gas densities 
become very low in the inner disk at late times), 
equation (7) shows that the fluffy grains will un¬ 
dergo slower radial drift, being more tightly cou¬ 
pled to the radial motion of the gas. However, 
even for fluffy grains, the dust migration is signif¬ 
icantly faster than that of the gas for radii > 10 
AU. This effect is much more marked for the com¬ 
pact grains, in which case the dust evolves orders 
of magnitude faster than the gas. Figure 2 shows 
the evolution of the accretion rate of the gas onto 
the star, plotted against the millimeter flux of the 
thermal dust (calculated for face-on disks at a dis¬ 
tance of 140 pc). The dashed lines demonstrate 
that this rapid evolution of compact grains vio¬ 
lates observational constraint (a) in §1, regard¬ 
ing the high incidence of millimeter emission in 
CTTSs. The dust disks of compact grains evolve 
to very low millimeter fluxes whilst the gas accre¬ 
tion rate is still high, thus generating a popula¬ 
tion of CTTSs without millimeter emission. The 
slower migration of the fluffy grains however keeps 
the millimeter flux above the detection limit (> 10 
mjy) until the gas accretion rate drops by a factor 
3 — 10. The variations in the disk parameters such 
as the initial disk mass (the initial surface density 
of the gas disk) and the viscosity result in rea¬ 
sonable agreement with the observed distribution 
of CTTSs in the plane of millimeter flux versus 
accretion rate. 

The dust evolution timescale is less than that 
of the gas for both types of grains (we would need 
to reduce p p by a further order of magnitude for 
the timescales to be comparable). This obviously 
means that constraint (b) in §1 is easily satisfied. 
Although the disk photoevaporation models pre¬ 
dict that gas resides at large radii (> 10 AU) for 
~ 10 7 yr after the cessation of accretion onto the 
star, such gas is, in these models, heavily depleted 
in dust. The predicted millimeter flux from the 
outer disk is orders of magnitude less than the 
current detection limit at the stage that the inner 
disk is cleared out. 

In summary, we find that models that combine 
viscous evolution, photoevaporation and gas-dust 
drag will readily satisfy constraints (b) and (c) in 
§1, but that constraint (a) is only satisfied in the 


case of fluffy grains (density p p < 0.1 g cm 3 ). 

3.2. Model Herbig AeBe Stars 

In this model, the central star is M = 2.5 Mg 
and L = 30 Lg (a typical value for Herbig AeBe 
stars; van den Ancker et al. 1998). The grav¬ 
itational radius, r g , beyond which ionized gas is 
unbound, increases by a factor 2.5. The initial 
value of the gas surface density is the same as the 
T Tauri star model, but the disk outer radius is 
increased to r ou t = 250 AU. The total disk mass is 
6.25 x 10 -2 Mg. Most significant, however, is the 
increase in ionizing luminosity. The observation of 
FUV (his < 13.6 eV) from young stars by Valenti 
et al. (2003) showed that continuum fluxes at 
A = 2257 A from Herbig AeBe stars are about two 
orders of magnitude larger than those of CTTSs 
on average. A stellar wind model by Bouret & 
Catala (1998) suggests that the ionizing fluxes of 
Herbig AeBe stars are of order of $41 = 10 2 — 10 4 . 
In our model, we assume 4*44 = 10 3 . 

The net photoevaporation mass loss rate, which 
is proportional to (Hollenbach et al. 

1994), is increased by a factor -\/2500 with re¬ 
spect to the T Tauri star model. The model thus 
evolves to the point that the inner disk is cleared 
out within a couple of Myr, an order of magnitude 
faster than the T Tauri star case. Consequently, 
the disk is much less depleted in the dust at the 
point that the inner disk is drained away. Figure 
3 demonstrates there is a substantial column in 
the dust in the outer disk at this stage. Evidently, 
once the wind truncates the disk at its inner edge 
(initially close to r 9 « 17 AU), the dust can no 
longer migrate toward the star and instead ac¬ 
cumulates in a dust ring close to the inner edge 
of the gas disk. We note that as the wind pro¬ 
gressively eats away at the outer disk, the inner 
edge of the gas disk grows to radii r r g (as in 
the models of Clarke et al. 2001) and that this 
is accompanied by the dust distribution becoming 
increasingly annular, as the dust is concentrated 
by inward migration, near the inner edge of the 
gas disk. The steep increase in the gas density 
and pressure with r near the inner edge induces a 
super-Keplerian velocity of the gas, and makes the 
inward velocity of the dust reverse in this region. 
This phenomena of dust accumulation is similar 
to those discussed by Klahr & Lin (2001) and by 
Haghighipour & Boss (2003a,b). It is intriguing 
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to note the similarity of this progression to the 
observed dust distributions in Herbig AeBe stars, 
such as the disk with a small inner hole in HD 
100546 (~ 10 AU; Bouwman et al. 2003) and the 
dust ring at a larger radius in HD 141569A (the 
inner edge ~ 175 AU; Weinberger et al. 1999; 
Augereau et al. 1999a; Clampin et al. 2003). 

We thus find that through changing the pa¬ 
rameters of the central star, but leaving the grain 
model unchanged, we can simultaneously account 
for the absence of long lived millimeter emission 
in WTTSs [constraint ( b ) in §1] whilst predicting 
outer dust disks of Herbig AeBe stars that have 
ceased to accrete. The essential difference is that 
the higher EUV flux of Herbig AeBe stars is able 
to significantly affect the disk on a much shorter 
timescale, so that the dust has not had time to 
migrate onto the star at the point that the inner 
disk becomes cleared. 

3.3. Other Parameters 

We now briefly consider how the above results 
depend on model parameters, while still retaining 
the assumption of a single sized grain population 
that is subject to neither grain creation nor de¬ 
struction processes. We find that we can vary the 
lifetime of the disk through varying a number of 
parameters (for example, the disk initial mass, ra¬ 
dius, or a viscosity parameter). The properties of 
the grains enter the evolution equations only via 
the combination p p s , which determines the degree 
of coupling between the grains and the gas (i.e., 
the stopping time in eq. [9]). We might thus ex¬ 
pect that the ratio of the evolution timescales of 
the dust and the gas should depend on this prod¬ 
uct, and therefore define a parameter 



This expectation is confirmed by Figure 4, 
which plots the outer dust disk lifetime against 
the inner gas disk lifetime for a range of disk 
models. The standard parameters of the model 
are the ones described in §3.1 for T Tauri disks, 
and the disk mass, the outer radius, the viscosity, 
and the dust parameter A are varied. Following 
Armitage et al. (2003), we use the gas accretion 
rate to determine the inner gas disk lifetime. The 
presence of gas accretion is generally considered 


as an observable indicator of having a gas disk. In 
our model, the gas lifetime is defined as the time at 
which the gas accretion onto the star drops to less 
than 10~ 10 Mq yr -1 . As shown by Clarke et al. 
(2001), when the inner disk is cleared out, gas ac¬ 
cretion vanishes. Hence, the gas lifetime defined 
here corresponds to the inner disk lifetime, al¬ 
though the outer disk might remain much longer. 
On the other hand, the outer dust disk lifetime 
is defined by the luminosity of 1.3 mm contin¬ 
uum. When the luminosity AttD 2 i'F 1/ drops below 
10 28 erg s -1 , the dust disk is defined to be re¬ 
moved. This criterion approximately corresponds 
to the detection limit of 2.5 mJy at D = 140 pc, 
which is the limit of the observation by Duvert et 
al. (2000). 

Clearly, the results lie in three bands of dust 
lifetime/gas lifetime, differentiated by the value of 
A. We thus see that, for the single sized grains of 
~ 1 mm, although the absolute lifetimes are depen¬ 
dent on the disk parameters such as the viscosity, 
the initial mass, and the radius, the relative life¬ 
times are just functions of A. Thus the fractional 
incidence of millimeter emission amongst CTTSs 
is controlled by A. This means the conclusion 
that millimeter sized grains must be rather fluffy 
(Pp ~ 0-1 § cm -3 ) in order to satisfy constraint 
(a) in §1 is a rather general one. 

3.4. Changing the Grain Size 

So far we have considered the case of grains 
of size 1 mm that are subject to neither replenish¬ 
ment nor destruction by coagulation or shattering. 
In CTTSs, the small values of opacity index at 
millimeter wavelengths, f3 < 1, (where oc A 3 ), 
imply that the grain size distribution extends to at 
least millimeter scales (Miyake & Nakagawa 1993; 
see however, e.g., Beckwith et al. 2000 for other 
interpretations). Thus, we do not further consider 
the case that the size s <C 1 mm. On the other 
hand, if the grain size is too large (s 1 mm), 
the opacity at millimeter wavelengths is reduced 
considerably (Miyake & Nakagawa 1993). 

If we now instead consider the case that s ~ 1 
cm, we have a problem, as discussed in detail by 
TL05. These grains undergo 10 times faster radial 
migration than millimeter sized grains and can¬ 
not survive over 10 6 yr. If the opacity at millime¬ 
ter wavelengths relies on the presence of centime¬ 
ter sized grains, then, in the absence of replen- 
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ishment, the timescale over which CTTSs would 
exhibit millimeter emission would be very short 
(less than 10 6 yr), in violation of constraint (a). 
Thus, in the absence of grain replenishment, we 
argue against the case s ~ 1 cm. 

4. Discussion 

4.1. Dust Grains Loss by the Photoevap¬ 
orating Wind 

The photoevaporating wind may carry the dust 
grains away from the disk. The gas drag force 
of the wind on a grain is F « ±nnom.HS 2 cf = 
rridnomiic 2 /(p p s), where is the spherical 
grain’s mass. For the Herbig AeBe star case, 
where M = 2.5 M 0 , r g = 2.5xl0 14 cm, $41 = 10 3 , 
and p p = 0.1 g cm -3 , the drag force is weaker than 
the gravity GMnid/r g , as long as s > 20 ^m. 
Hence, if the dust grains have grown larger than 
20 pm, most of the dust mass remains in the disk 
even after the gas disk has completely evaporated. 
In this case, it is expected that a ring-like dust disk 
remains after the gas dispersal, like the HR 4796A 
dust-ring (Schneider et al. 1999). Further evolu¬ 
tion of the dust ring during or after the gas dis¬ 
persal likely alters its morphology, for example, by 
radiation pressure (Klahr & Lin 2001; Takeuchi & 
Artymowicz 2001), by collisions of grains (Kenyon 
& Bromley 2004), and by unseen planets (Wyatt 
et al. 1999). (See also Augereau et al. 1999b; Li 
& Lunine 2003; Currie et al. 2003; Chen & Kamp 
2004, for modeling of the HR 4796A disk based 
on its emission.) Detailed discussion on such later 
evolution stages are, however, beyond the scope 
of this paper. 

For T Tauri disks of <l> 4 i = 1, the drag force 
by the wind is much weaker, and thus even 1 pm 
grains can resist the wind. Hence, grain loss by 
the inward migration is much more important in 
T Tauri disks. 

4.2. The Possible Role of Grain Replen¬ 
ishment 

We concluded in §3.3 and §3.4 that the bulk of 
the opacity at millimeter wavelengths is provided 
by grains of sizes similar to or less than about 1 
mm in order to satisfy constraint (a). Regarding 
the fraction of CTTSs with detectable millime¬ 
ter emission, we need to retain a stock of small 
(s < 1 mm) grains in the disk over several Myr. 


We have shown in §3 that such grains could avoid 
excessively rapid migration if they were rather 
fluffy, but there is also a possibility that they are 
otherwise depleted by coagulation to larger grain 
sizes. This possibility was considered in TL05. In 
this subsection, we note that this problem would 
be circumvented if the stock of grains of s < 1 
mm were replenished, either by continued infall of 
small grains from a circumstcllar envelope or by 
shattering of larger bodies. 

We cannot rule out that this is the case. How¬ 
ever, we note that although this would certainly 
satisfy constraint (a), it would, without fine tun¬ 
ing, fail to satisfy constraint ( b ). In other words, if 
opacity at millimeter wavelengths were sustained 
in this way throughout the CTTS phase, it would 
over produce the millimeter emission from the 
outer disk in the WTTS phase. The only way that 
this could be avoided would be if, coincidentally, 
the transition from the CTTS phase to the WTTS 
phase also turned off the replenishment process. It 
is not clear, however, why this should be the case. 

5. Conclusions 

We have shown that the disk model that com¬ 
bines viscous evolution, photoevaporation and the 
differential radial motion of grains and gas can ac¬ 
count for the three observational constraints (a)- 
(c) set out in §1. We particularly highlight the fact 
that the same model of the grains can simultane¬ 
ously account for the fact that the non-accreting 
T Tauri stars (WTTSs) do not possess detectable 
millimeter emission, whereas some of their non¬ 
accreting counterparts in the Herbig AeBe mass 
range show prominent dust disks at large radii. 

The difference between the two types of behav¬ 
ior is due to the fact that Herbig AeBe stars suf¬ 
fer much stronger photoevaporative mass loss from 
the disk (due to the stronger EUV flux) and this 
causes accretion onto the star to cease at a much 
younger age. This relatively rapid draining of the 
inner gas disk means that there is insufficient time 
for the grains to have migrated onto the star. Con¬ 
sequently, the grains are stranded in the outer disk 
(i.e., at r > r g ). As the outer disk is slowly eroded 
by photoevaporative mass loss, the inner hole in 
the disk grows, while at the same time the dust be¬ 
comes increasingly concentrated toward the inner 
edge of the residual gas disk. Thus the appearance 
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of dust rings at large radii (tens of AU), which 
is sometimes attributed to the action of planets, 
is a natural consequence of photoevaporation and 
grain migration. In T Tauri stars, by contrast, 
the photoevaporative flow is much gentler. By 
the time this flow has created an inner hole in 
the gas disk, the dust has all migrated onto the 
star, and so WTTSs are not generally millimeter 
wavelength sources. It is intriguing in this regard 
that the small number of T Tauri stars with op¬ 
tically thin inner disks which yet show millimeter 
emission from the disk at large radii are objects 
with unusually high ultraviolet fluxes (Bergin et 
al. 2004). This supports our conjecture that the 
timescale for mass loss due to photoevaporation, 
relative to the timescale for grain migration, is an 
important factor in the formation of these systems. 

We find that the hardest of the observational 
constraints to satisfy is that the millimeter flux in 
T Tauri stars should persist for a good fraction 
of the time that the star spends as a CTTS. This 
places a constraint on the relative timescales for 
dust and gas evolution, which, in our models, is 
determined almost entirely by the product of the 
grain density and the grain size, p p s, since this 
controls the degree to which the grain motions are 
coupled to the gas. We find that in order for the 
millimeter flux to persist long enough we need to 
invoke a population of grains that is rather fluffy 
( P P < 0.1 g cm -3 ). 

Finally, we note a clear observational prediction 
of the type of model that we have described here. 
In these models, WTTSs are systems where the 
inner disk has drained onto the star, having been 
starved of resupply by photoevaporative mass loss 
from the outer disk. WTTSs then retain disk gas 
at large radii for considerable periods (Clarke et al. 
2001), but we argue here that the dust content of 
such disks should be low, due to the prior radial 
migration of the dust. In this way, we circum¬ 
vent the problem that WTTSs show no detectable 
dust emission in general. We however predict that 
these systems should retain a significant column 
density in gas during the WTTS phase (see Fig. 
1). We point out that detection of gas diagnostics 
in WTTSs at large radii would provide striking 
confirmation of this type of model. 
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Fig. 1. — Evolution of the gas and dust surface 
densities of a T Tauri disk. The dust surface 
densities are multiplied by the gas-to-dust ratio 
(/dust = 100) and shown by the solid lines. The 
dashed lines show the gas surface densities, (a) 
Compact grains, p p = 1 g cm -3 . Most of the 
dust grains have disappeared by 10 6 yr. (b) Fluffy 
grains, p p = 0.1 g cm -3 . The dust grains remain 
over 10 6 yr, but they cannot survive until photo¬ 
evaporation begins to clear the inner gas disk at 
1.6 x 10 7 yr. 


x x Ivl =0.25 M 



Fig. 2.— Gas accretion rate vs. dust thermal 
flux at a wavelength of 1.3 mm. The solid lines 
are for disks with fluffy grains ( p p = 0.1 g cm -3 ), 
while the dashed lines are for compact grains 
(p p = 1 g cm -3 ). The evolution of three mod¬ 
els [the standard model, the massive disk model 
of M g = 0.25 Mq (10 fold increase in the ini¬ 
tial gas density), and the high viscosity model of 
a = 10 -2 ] is plotted. The 1.3 mm fluxes are calcu¬ 
lated for face-on disks at a distance of D = 140 pc. 
Crosses are the observed values of T Tauri stars. 
The 1.3 mm fluxes are taken from Beckwith et al. 
(1990) and Osterloh & Beckwith (1995), and the 
gas accretion rates are taken from Hartmann et al. 
(1998). Squares show the upper limits of 1.3 mm 
fluxes. 
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Fig. 3.— Same as Fig. 1, but for a Herbig AeBe 
disk. The grains arc fluffy ( p v = 0.1 g cm -3 ), and 
the ionizing photon luminosity is $41 = 10 3 . 
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Fig. 4.— Lifetimes of the gas and dust disks of T 
Tauri star models. The upper group of the thin 
symbols are for A = 10~ 2 , the thick symbols are 
for A = 0.1, and the lower thin symbols are for 
A = 1. Asterisks (*) are for the standard gas disk 
with M g = 2.5 x 10 -2 M 0 , r out = 100 AU, and 
a = 10 -3 . Crosses (+) are for disk masses, M g = 
2.5 x 10" 3 , 8 x 10" 3 , 8 x 10" 2 , and 2.5 x 10" 1 M Q 
from the left. Triangles (A) are for disk radii, 
Aut = 30 AU and 300 AU from the left. Diamonds 
(O) are for viscosities, a = 10 -2 and 10 -4 from the 
left. The dotted line represents that the lifetimes 
of the gas and the dust are the same. 
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